Extrinsic Entwined with Intrinsic Spin Hall Effect in Disordered Mesoscopic Bars 
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We show that pure spin Hall current, flowing out of a four-terminal phase-coherent two- 
dimensional electron gas (2DEG) within inversion asymmetric semiconductor heterostructure, con- 
tains contributions from both the extrinsic mechanisms (spin-orbit dependent scattering off impuri- 
ties) and the intrinsic ones (due to the Rashba coupling). While the extrinsic contribution vanishes 
in the weakly and strongly disordered limits, and the intrinsic one dominates in the quasiballistic 
limit, in the crossover transport regime the spin Hall conductance, exhibiting sample-to-sample large 
fluctuations and sign change, is not simply reducible to either of the two mechanisms, which can be 
relevant for interpretation of experiments on dirty 2DEGs [V. Sih et at, Nature Phys. 1, 31 (2005)]. 

PACS numbers: 72.25.Dc,71.70.Ej, 73.23.-b 
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Introduction. — Recent experimental discovery of the 
spin Hall effect in three-dimensional n-type semiconduc- 
tors and two-dimensional hole 01 or electron gases ^3] 
provides deep insight into the relativistic effects in solids 
which manifest through spin-orbit (SO) couplings 0|, as 
well as into the nature of spin accumulation and spin 
fluxes in systems whose physics is governed by the SO in- 
teractions. The flow of conventional unpolarized charge 
current in the longitudinal direction through such sys- 
tems leads to spin flux in the transverse direction which 
deposits nonequilibrium spin accumulation on the lateral 
sample edges, as observed in experiments 0,0, 0I- More- 
over, the anticipated experiments should demonstrate 
flow of spin into the transverse electrodes attached at 
those edges 0, thereby offering a semiconductor analog 
of the Stern-Gerlach device as a source of spin currents 
or flying spin qubits where spatial separation of spin-f 
and spin-j electrons does not require any cumbersome 
external magnetic flelds. 

However, the effect observed in electron systems P, Q 
is rather small, and its magnitude cannot be tuned easily 
since it is determined by material parameters which gov- 
ern the SO dependent scattering off impurities 0, Q, 0] 
deflecting a beam of spin-| (spin-j) electrons predomi- 
nantly to the right (left) . In contrast to such extrinsically 
(i.e., impurity) induced spin Hall effect, recent theories 
have argued for several orders of magnitude larger spin 
Hall currents due to the intrinsic SO couplings capable 
of spin-splitting the energy bands and inducing trans- 
verse SO forces |^ in bulk semiconductors in the clean 
limit |9|,|l3, as well as in ballistic mesoscopic multitermi- 
nal nanostructures made of such materials |^ . Such 
mechanisms were invoked to explain |l2| observation of 
two orders of rnagnitude larger spin Hall accumulation in 
2D hole gases |3, whose ma grii tude can still be affected 
by the effects of the disorder For example, in infinite 
2DEGs with linear in momentum SO coupling arbitrar- 
ily small amount of impurities completely suppresses the 
spin Hall effect while increasing disorder only grad- 
ually diminishes the signatures of the spin Hall effect in 
mesoscopic 2DEG structures ^3, . 



Despite these advances, the intrinsic vs. extrinsic 
debate on the origin of the detected signatures of the 
spin Hall effect persists 0| , closely mimicking sev- 
eral decades old discourse on the Karplus-Luttinger in- 
trinsic vs. skew-scattering -I- side-jump (i.e., shift of 
the scattered wave packet) possible explanations of 
the anomalous Hall effect observed in ferromagnetic sys- 
tems ,19] . Surprisingly enough, no theory has been con- 
structed to describe spin current in the transition from 
quasiballistic to dirty transport regime in multiterminal 
Hall bars employed in experiments where, in principle, 
both extrinsic and intrinsic mechanisms could be equally 
important. Also, at low temperatures, quantum correc- 
tions to the semiclassical description 0,0 of the extrinsic 
spin Hall effect, such as the mesoscopic fluctuations and 
localization effects due to the interplay of randomness 
and quantum coherence, can play an essential role for in- 
terpreting experiments on meso- and nano-scale samples. 

The demand for such theory is further provoked by 
very recent spin Hall experiment on dirty two-terminal 
2DEG within AlGaAs quantum well which accommo- 
dates both the SO scattering off impurities and a small 
Rashba SO coupling due to weak structural inversion 
asymmetry. The latter could give rise to an intrinsic 
contribution of a similar magnitude as the extrinsic one 
for not too small ratio Aso^/fi 0,0] (Aso is the spin- 
splitting of quasiparticle energies and fi/r is the disorder 
induced broadening of energy levels due to flnite trans- 
port scattering time r) — the parameters measured in 
Ref. 3| place its Hall bar into the regime Aso'^'/fi- ^ 10^^. 

Here we construct a quantum transport theory of the 
spin Hall fluxes driven by a generic Haniiltonian 
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which encompasses both the intrinsic (third term) and 
the extrinsic (fourth term) SO coupling-induced ef- 
fects. Within this framework, we obtain in Fig. Q] 
the spin Hall conductance in the transverse ideal (i.e., 
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FIG. 1: (Color online) The disorder-averaged typical spin 
Hall conductance of a four-terminal phase-coherent 2DEG 
nanostructure of the size 100a x 100a (a ~ 2 nm), which 
is governed by both the SO scattering off impurities (setting 
the transport scattering time r) of strength Aso and the in- 
trinsic Rashba SO coupling responsible for the spin-splitting 
Aso = Aatsokp- Following measured parameters of 2DEG 
in Ref. 0, we use Aso = 0.005, tso = for the "extrinsic"; 
Aso ~ 0, tso = O.OOSto, for the "intrinsic"; and Aso ~ 0.005, 
tso = O.OOSfo for the "extrinsic -f- intrinsic" spin Hall con- 
ductance GaH ~ li/iVi — V4) of the pure spin current /| in 
the transverse ideal lead 2. 



spin and charge interaction-free) leads of a four-terminal 
quantum-coherent Hall bar containing such 2DEG. While 
the intrinsic SO coupling term has general Zeeman form 
—gUBO' ■ B(p)/2, its effective magnetic field B(p) is mo- 
mentum dependent and does not break the time-reversal 
invariance. For the Rashba coupling, which arises due 
to structural inversion asymmetry (of the 2DEG confin- 
ing electric potential along the z-axis and differing band 
discontinuities at the heterostructure quantum well in- 
terface I3j), B(p) — ~{2a/ g^B){p x z) where z is the 
unit vector orthogonal to 2DEG. The corresponding spin- 
splitting of the energy bands at the Fermi level (hkp is 
the Fermi momentum) is Ago = (7/^b|B(p)| = 2akp. 

The Thomas term A (it x p) • VVdis(i") is a relativis- 
tic correction to the Pauli equation for spin-^ particle 
where A/fi = -fi^/Amlc^ « -3.7 x 10^^ (c is velocity 
of hght) for electron in vacuum with mass mo. While 
this vacuum value would lead to hardly observable spin 
Hall effect, in solids A can be enormously renormalized 
by the band structure effects due to the strong crystal 
potential — for example, X/h = 5.3 in GaAs 0] gen- 
erates the extrinsic spin Hall conductivity af^"^ = jy / 
whose recently computed value 0, could account for 
the magnitude of the observed spin Hall accumulation in 
dirty 3D electron systems of Ref. P|. In the inversion 
symmetric systems B(p) = 0, so that pure spin current 
jy (carrying z-axis polarized spins along the y-axis) is 
driven solely by the Thomas term with A 7^ in the 
presence of longitudinal external electric field Ex ■ 



FIG. 2: (Color online) The full distribution function P{Gsh) 
of the spin Hall conductance Gsh fluctuating from sample 
to sample in an ensemble of 10000 impurity configurations 
within 2DEG with the extrinsic SO scattering Aso = 0.005 
and the intrinsic Rashba SO coupling tso = O.OOSto. The in- 
set shows the corresponding distribution function of voltages 
V2 which have to be applied to the top transverse lead (to- 
gether with voltage V3 on the bottom lead) to ensure that no 
charge current I2 = I3 = accompanies pure spin current /| . 



and the side-jump a^fj > contributions, is a singu- 
lar function as the disorder strength is decreased to zero 
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The extrinsic spin Hall conductivity cr^^"^ 
which is the sum ^, ^ of the skew-scattering 
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T while fjf]^ is r-independent 0, H, E3- 
Therefore, in the clean limit this divergence is expected 
to be cut off at r ~ H/Ago, where Aso is due to 
the k'^ Dresselhaus term in the absence of Rashba cou- 
pling 0- In contrast to this unphysical behavior of 
bulk transport quantities. Fig. ^ shows that the extrin- 
sic spin Hall conductance of the four-probe mcsoscopic 
bar, G°^'' = GsniX / 0,a = 0) goes to zero GsH 
in the clean limit Ago S> h/r, as well as in the very 
{kpi ^ 1) dirty limit Aso fi/T where semiclassical 
theories break down. The full spin conductance is de- 
fined as GsH = li/iVi ~ V4) for transverse spin current 
^2 = ^i^2~^2) of z-polarized spins driven by the voltage 
bias Vi —Vi [see inset of Fig. p. 

Although we find that the disorder-averaged extrinsic 
spin Hall conductance (Gf^^) = {GsHiX,a = 0)) de- 
pends linearly on A, as is the case of the extrinsic bulk 
spin Hall conductivity cr®^'' cx A, Fig. [21 demonstrates 
that in phase-coherent nanostructures one also has to 
take into account violent mesoscopic fluctuations on the 
top of semiclassical background — they affect the change 
of the sign of GsH from sample to sample and generate 
wide distribution of conductances P{Gsh)- Thus, we use 
the typical value {Gsh) = (G^ff) typical, at which the 
corresponding full distribution function P{Gsh) reaches 
its maximum (see Fig. |5] for illustration), to characterize 
an ensemble of mesoscopic spin Hall bars differing from 
each other in the impurity configuration. Note that both 
(Gs//)typicai and (Gsfl-)mean are positivc within the dirty 
regime of recent experiments meaning that the spin 

Hall current /| flows from the top to the bottom trans- 
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verse lead because of spin-t being deflected to the right, 
in accord with reported sign of lateral spin accumulation 
in two-terminal microstructures P, Q . 

Extrinsic + intrinsic total spin Hall currents in ter- 
minals of mesoscopic bars. — Theoretical analysis of the 
spin Hall effect in bulk systems yields asH by comput- 
ing (via the Kubo formula d, 0, 0, Q| or kinetic equa- 
tions |2iiJ) the pure spin current density jy, as the ex- 
pectation value of some plausibly introduced spin cur- 
rent operator jy , in linear response to external elec- 
tric field Ex penetrating an infinite semiconductor. This 
procedure is well-defined in inversion symmetric systems 
(such as A 7^ 0, a = 0) where SO contribution to the 
impurity potential decays fast and spin current is a con- 
served quantity 0]. However, non-conservation of spin 
in spin-split systems with B(p) ^ leads to ambiguity 
in the definition of spin current and the corresponding 
o'sH E3- That is, the spin density relaxes not only due 
to the spin flux but also due to spin precession, thereby 
requiring careful treatment of boundary effects on the 
observed spin accumulation ^ • 

To generate analogous phenomenology in multiter- 
minal Hall bars we have to employ at least four ex- 
ternal leads where passing unpolarized charge current 
through the longitudinal electrodesgenerates spin cur- 
rent in the transverse electrodes |M 0|. The charge 
currents in the terminals are described by the same 
multiprobe Landauer-Biittiker scattering formulas em- 
ployed for mesoscopic quantum Hall bars ^ 
J2q Gpq{Vp — Vq), which relate current in probe p to volt- 
ages Vq in all probes attached to the sample via conduc- 
tance coefficients Gpq that intrinsically take into account 
all interfaces and boundaries. To ensure that transverse 
spin Hall current is pure, one has to apply transverse 
voltages V2 and V3 obtained from these equations for 
h ~ I3 = 0. The extension of the scattering theory 
of quantum transport to total spin currents in the 
terminals 1^ — j^i^p — Ip), which are always genuine 



nonequilibrium response (in contrast to spin current den- 
sity which can be non-zero even in equilibrium ]&) and 
conserved in the ideal [B(p) = 0] leads, yields 
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where G% - + G],^ - G^^J - Gj,^, G™' - G" + 
Gql — GJ^ — G^^, and the standard charge conductance 
coefficients are Gpq = G^J -\- Gl\ -t- G^J -h G^\. The 
spin-resolved conductances are obtained from the trans- 
mission matrices t^'' between the leads p and q through 
the Landauer-type formula G^^' = x E»,i=i Itfj^tr' 
where J2i j=i l^^jaa'l'^ ^he probability for spin-a' elec- 
tron incident in lead q to be transmitted to lead p as 
spin-(T electron and i,j label the transverse propagating 
modes in the leads. Solving Eq. for /| in the top 
transverse lead of the Hall bar in Fig. ^ yields a general 
formula 11 ll 
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where we choose the reference potential ¥4 = and the 
2;-axis is selected as the spin quantization axis for ti i- 

The consequences of the scattering approach-based 
four-terminal spin Hall conductance formula Eq. lO can 
be worked out either by analytical means, such as the 
random matrix theory applicable to "black-box" disor- 
dered or chaotic ballistic structures, or by switching to 
nonperturbative real^spin space Green functions 5] to 
take into account microscopic Hamiltonian which models 
the details of scattering and SO coupling effects in finite- 
size samples attached to many external probes. For this 
purpose, we represent the Rashba 111 and the Thomas 
term (2^ in Eq. JQ) in the local orbital basis 



m,CT (mm') crcr' m,a/3 ij 



The first term models isotropic short-range spin- 
independent static impurity potential where Em G 
[—W/2,W/2\ is a uniform random variable. The sec- 
ond, Rashba term, is a tight-binding type of Hamiltonian 
whose nearest-neighbor (mm') hopping is a non-trivial 
2x2 Hermitian matrix t^i^ — (tm^n')^ ^"^^ spin space 

_{ -tola - itso^y (m = m' -f e^;) , , 

\ -toIs+iiso-T. (m = m' + ej,). ^''> 



The strength of the SO coupling is measured by the pa- 
rameter tso — a /2a {Is is the unit 2x2 matrix in the spin 
space) which determines the spin-splitting of the band 
structure Aso = lafso^F- A direct correspondence be- 
tween the continuous effective Hamiltonian Eq. and 
its lattice version Eq. 10} is established by selecting the 
Fermi energy {Ep = — 3.5to in the rest of the paper) 
of the injected electrons to be close to the bottom of 
the band where tight-binding dispersion reduces to the 
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parabolic one, and by using to = It? / {2'm* a?) for the or- 
bital hopping which yields the effective mass to* in the 
continuum limit. The labels in the third (Thomas) term, 
which involves second neighbor hoppings, are the dimen- 
sionless extrinsic SO scattering strength Aso = Afi./ (4a^), 
eijz stands for the Levi-Civita totally antisymmetric ten- 
sor with i,j denoting the in-plane coordinate axes, and 
zy, 7 are the dummy indices taking values ±1. For 2DEG 
of Ref. 0] these SO coupling parameters take values 
tso — O.OOSto and Aso — 0.005, assuming conduc- 
tion bandwidth ~ 1 eV and using the effective mass 
TO* = 0.074TOO 0. 

Mesoscopic fluctuations and scaling of the spin Hall 
conductance. — Akin to the dichotomy in the description 
of the quantum Hall effect [2^, the characterization of 
the spin Hall effect in terms of cr.sff is (assuming "reli- 
able" definition of spin current [l7|) applicable at suf- 
ficiently high temperatures ensuring the classical regime 
where dephasing length is much smaller than the sam- 
ple size and conductivity can be treated as a local quan- 
tity. On the other hand, in mesoscopic samples of size 
L < Lrj, nonlocal quantum corrections to semiclassical 
background and presence of external measuring circuit 
have to be taken into account through GsH- 

For example, quantum coherence and stochasticity in- 
troduced by the impurities will lead to sample-to-sample 
fluctuations of GsH , so that lack of self-averaging requires 
to study its full distribution function P{Gsh), as demon- 
strated by Fig. 12 Nevertheless, the main part of the 
distribution is best characterized by the typical (i.e., geo- 
metric mean) conductance in Fig. 13 rather than the 
usual arithmetic mean value affected by the long tails 
of P{Gsh)- The inset of Fig. |21 also reveals sample-to- 
sample fluctuations of the voltages that have to be ap- 
plied to transverse leads to ensure pure spin Hall current 
I2' — ^I^' (note that extrinsic SO scattering generates 
only the z-polarized spin Hall current studied here, while 
Rashba SO coupling is also responsible for non-zero spin 
currents with in-plane polarized spins with no net 
charge current I2 = I3 = in the transverse direction. 
In perfectly clean and symmetric bars this condition is 
satisfied by applying V2 = V3 ^ {Vi - V^)I2 

Although the samples in Fig. El showing the scaling of 
{Gsh){L), are of the size L <1 fim (compared to 100 ^m 
sizes of samples required for optical scanning of the spin 
Hall accumulation 0, Q), their properties can be used 
to account for spin Hall quantities of microstructures 
viewed as a classical network of resistors, each of the size 
X L^. The conductance {Gsh){L — L^) of each resis- 
tor can be read from Fig. |21 as computed in fully phase- 
coherent regime with localization length ^ ^ L > Lip. 
Even though inelastic scattering will destroy quantum- 
coherent nature of transport between resistors, thereby 
pushing it into the semiclassical regime 0, Q| , the spin 
Hall conductance of the whole network (i.e., of a macro- 
scopic 2DEG bar at high temperatures) is still equal to 




2DEG Size L (a) 

FIG. 3: (Color online) Scaling of the typical spin Hall con- 
ductance with the size of the 2DEG characterized by the ex- 
trinsic SO scattering strength Aso = 0.005 and the intrinsic 
Rashba SO coupling tso = O.OOSio which, together with the 
disorder setting the mean free path £ — vft = 34a, deter- 
mine Asor/h = 0.2. For comparison, the "extrinsic" curve 
corresponds to tso = = Aso. 



{Gsh){L — Lp) in two-dimensions. 

Conclusions. — Our principal result in Fig. ^ demon- 
strates that the spin Hall effect in four-terminal disor- 
dered 2DEGs can be affected by both the SO scattering 
of impurities and the Rashba SO coupling inducing the 
energy spin-splitting Aso- The former dominates GsH — 
G^^'' in the dirty regime AgoT/h < 10~^, while the latter 
accounts for GsH — G^^^'^ in the clean regime Asor/ft > 1 
where it could be tuned via the gate electrode [2^ to two 
orders of magnitude larger value (for presently achiev- 
able Rashba couphng strengths tgo ~ 0.01 than the 
maximum value of G"^' reached in the dirty limit. In the 
crossover regime 10~^ < Asot/?i < 1, our nonpertur- 
bative quantum transport analysis unveils an interplay 
of the extrinsic and intrinsic mechanisms responsible for 
GsH ^ Gf^^ + G^J'^^. These predictions offer a clear guid- 
ance toward nanofabrication of high-mobility devices in 
order to increase the magnitude of the spin Hall current, 
as well as to evade large sample-to-sample fluctuations 
of GsH and its sign change in Fig. |21 characterizing the 
low-temperature quantum spin Hall transport in dirty 
phase-coherent bars. 
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